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We point out that Bose enhancement in a hadronic wave function generically leads to correlations 
between produced particles. We show explicitly, by calculating the projectile density matrix in the 
Color Glass Condensate approach to high-energy hadronic collisions, that the Bose enhancement 
of gluons in the projectile leads to azimuthal collimation of long range rapidity correlations of the 
produced particles, the so-called ridge correlations. 


I. INTRODUCTION 

The ridge structure observed in high multiplicity p-p 
[l] and p -Pb Q collisions at the Large Hadron Collider 
(LHC) triggered an intense activity aimed at understand¬ 
ing the possible physical origin of correlations between 
emitted particles. Two basic ideas have been put for¬ 
ward in this context (see others in Q). 

According to one idea, the origin of the correlations 
is the same as in similar ridge correlations observed ear¬ 
lier in heavy ion collisions at the Relativistic Heavy Ion 
Collider fj and the LHC Q. Namely, the angular colli¬ 
mation is due to flow effects in the final state The 
qualitative features of the high multiplicity p-p and p -Pb 
data, including the dependence on masses of produced 
particles, are well described by the hydrodynamic-based 
models. It is nevertheless challenging to explain how the 
spatially small system produced in the final state in p- 
p collisions can sustain the collective behavior necessary 
for local equilibration. 

The second suggestion is that the final state corre¬ 
lations carry the imprint of the partonic correlations 
that exist in the initial state. Three different variants 
of such initial state effects have been discussed in the 
literature: local anisotropy of target fields spatial 
variation of partonic density [U and finally the “glasma 
graph” contributions to particle production @ within the 
Color Glass Condensate (CGC) approach to high-energy 
hadronic scattering. While the physical origin of the 
first two effects is quite clear, the physics behind glasma 
graphs has not been elucidated in the literature. On the 
other hand, numerical calculations based on the glasma 
graph approach have been very successful in reproducing 
the systematics of ridge correlations 0 - It is therefore 
important to understand the physics that underlies these 
numerical results. 

The purpose of this Letter is to point out that there ex¬ 
ists a general quantum mechanical mechanism that leads 
to positive correlations of emitted particles with simi¬ 
lar quantum numbers. It is operative when the wave 


function of an incoming hadron is dominated by bosons 
(gluons), and is due to Bose enhancement in this wave 
function. In the next section, after recalling the basic 
derivation, we will show that this is precisely the physical 
mechanism that underlies the glasma graph calculation of 
hadron production in p-p and p -A collisions. The mecha¬ 
nism itself is however more general, has been widely used 
for identical mesons in heavy-ion collisions, see e.g. 0- 
Analogously, for fermions in the initial state one expects 
the opposite effect, namely Pauli blocking. In the final 
section, after discussing our results, we briefly address 
the question of which final state observables could be 
sensitive to the initial state Pauli blocking. 

To avoid confusion, we stress that by Bose enhance¬ 
ment we do not mean the Hanbury-Brown-Twiss (HBT) 
correlations between emitted particles, which arise due to 
emission from a large source comprised of many incoher¬ 
ent emitters. We rather mean the effect of Bose statistics 
that enhances the probability to find two identical bosons 
with the same transverse momentum in the incoming pro¬ 
jectile wave function before the collision. Although this 
initial state enhancement is for two incoming bosons with 
identical transverse momentum, the momenta of the two 
are modified differently by the interaction with the tar¬ 
get. Thus after scattering the two bosons emerge in the 
final state with different momenta. The correlation be¬ 
tween the directions of the two momenta is nevertheless 
preserved in some range of kinematics, see later. 


II. GLUON PRODUCTION AND BOSE 
ENHANCEMENT 

A. Basics of Bose enhancement 

The prototypical textbook calculation of Bose 
enhancement proceeds as follows 0. Consider a 
state with fixed occupation numbers of N species 
of bosons at different momenta, | {^(p)}) = 

n, P ^2 (p) , (4 (p)/vv) n '(p)\o), i = iwith 

a finite volume V and periodic boundary conditions so 
that momenta are discrete. The state is translationally 




invariant with mean particle density 

n ee ({n\p)}\a^(x)a i (x)\{n\p)}) = £ n\p). (1) 

i,p 

Hereafter we take ')2 [l « / d 3 p/( 2tt) 3 . The 2-particle 
correlator in coordinate space is 

D{x,y ) = ({n(p)}\a ] \x)a h (y)a\x)a° {y)\{n{p)}). (2) 


This is calculated by going to momentum space, where 
the operator averages are simple: 

({n{p)}W\p)a^{q)a\l) a >{m)\{n{p)}) (3) 

= ({n(p)}\S(p - l)S(q - m)a^{p)a l {p)a f:l {q)a J {q) 

+ S(p — m)S(q - l)a} l (p)a? (p)a h (q)a l (q)\{n(p)}) 

m S(p - l)6(q - m)Y^n l {p)'^2 n ' (?) 

* 3 

+ S(p - m)S(q - l)Y^n l {p)n\q), 

i 

where we have neglected the terms where all momenta 
are equal, which are suppressed by a phase space factor. 
Using this, the result for D{x , y) reads 


D(x,y) =n 2 + 

i 


d P Jp(. x ~iO nHr,) 
(2tt)3 6 U {P) 


( 4 ) 


The last term expresses the Bose enhancement. It van¬ 
ishes when the points are very far away, and gives 0(1/N) 
enhancement when the points coincide. The 0(1/N) sup¬ 
pression of the second term relative to the first one is due 
to the fact that the second term contains a single sum 
over the species index. The physics is that only bosons 
of the same species are correlated with each other. Tech¬ 
nically the origin of this additional contribution is the 
“wrong contraction” term in eq. 0- 
The Bose enhancement is a generic phenomenon, and 
is not tied specifically to the state with fixed number of 
particles. An overwhelming majority of pure states or 
quantum density matrices exhibit Bose enhancement at 
some degree. There is however one type of states that do 
not exhibit such behavior, notably classical-like coherent 
states. Consider a coherent state 


|6(x)) = exp{* 


d 3 x & < (z)(a < (») + a ti (a;))}|0>. 


( 5 ) 


A trivial calculation in this state gives 
(b(x)\a^(x)a l (x)\b(x)) = b l (x)b l (x), (6) 

(' b{x)\a ]l {x)a ]3 {y)a\x)a J {y)\b{x)) = b l (x)b l {x)V (y)b J (y ), 

so D(x,y) = n{x)n{y). Thus, in order to exhibit Bose 
enhancement, a state has to be nonclassical. 


B. Gluon production via glasma graphs 

As stated above, we want to demonstrate that the 
angular collimation arising from the glasma graph cal¬ 
culation owes its existence to the Bose enhancement in 


the projectile wave function. Following [IS in , we con¬ 
sider the calculation of inclusive two particle production 
and assume local parton-hadron duality, namely that at 
a given momentum the number of produced hadrons is 
proportional to the number of produced gluons. 

The graphs that contribute to this observable can be 
conveniently represented in terms of averages of gluon 
creation and annihilation operators in the incoming pro¬ 
jectile wave function (see supplementary material and 
jl3l|). They come in three varieties, see fig. [I] Type 
A graphs give the contribution whereby two gluons from 
the incoming projectile wave function scatter indepen¬ 
dently on the target. The incoming gluons have trans¬ 
verse momenta Aqand k 2 respectively. While propagating 
through the target the first particle picks up transverse 
momentum p — k\ and the second particle picks up trans¬ 
verse momentum q — k 2l so that the outgoing particles 
have momenta p and q. Type B and C graphs from the 
projectile point of view are “interference graphs”, in the 
sense that the final state gluon with momentum p comes 
from the projectile gluons with different momenta in the 
amplitude and complex conjugate amplitude. Type B 
and C diagrams contain leading contributions that can 
be reinterpreted as Type A but with gluons originating 
from the target rather than from the projectile, and ad¬ 
ditionally subleading contributions, including those that 
lead to HBT correlations Q- This is explained in detail 
in the Appendix. Therefore, in the following we will only 
discuss those of Type A, keeping in mind this comple¬ 
mentary interpretation of the leading pieces of Type B 
and C. 

The Type A contribution to double inclusive gluon pro¬ 
duction can be written as 

C f (in\al\k 1 )al\k 2 )a^(k 1 )a l b (k 2 )\in) (7) 

Jk \,&2 

h,i h.k 

x 

p2. 

where | in) is the wave function of the incoming projec¬ 
tile, C is a constant, N(p — k) is the probability that the 
incoming gluon with transverse momentum k acquires 
transverse momentum p after scattering and, hereafter, 
we use the notation f k = J ■ This scattering proba¬ 
bility is of course determined by the distribution of target 
fields (within the glasma graph calculation, the scatter¬ 
ing of the two gluons is independent). We have also as¬ 
sumed that the target wave function is translationally 
invariant, so that the momentum transfer is the same to 
the gluon in the amplitude and complex conjugate ampli¬ 
tude. The last assumption does not allow one to discuss 
the correlation mechanisms proposed in [3, H»|] within this 
framework. 

Also note that we have not indicated the rapidity vari¬ 
able on the gluon creation and annihilation operators. 
Within the glasma graph calculation the gluon produc¬ 
tion is rapidity independent. Rapidity dependence be¬ 
comes significant only when the rapidity difference be- 
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FIG. 1. Glasma graphs for two gluon inclusive production before averaging over the projectile color charge density p. Black 
blobs denote vertices which involve two momenta, e.g. W — p 1 k 3 1 /p 2 , and dashed lines the cuts. For details see the text. 


tween the observed hadrons becomes large, Ar/ ~ l/a s . 
The origin of this independence is that the CGC hadronic 
wave function is approximately boost-invariant. In fact, 
only the rapidity independent mode of the gluon field is 
large in the wave function of the fast hadron, and only the 
creation operators of this one rapidity mode are relevant 
to the discussion of correlations. 

Thus, the creation and annihilation operators enter¬ 
ing the above equation are the original gluon operators 
integrated over rapidity, 

ai(k) = -L[ p-ai(rt,k). ( 8 ) 

VY J\ri<Y/2\ 27r 

Here the rapidity interval Y/2 is arbitrary, but large 
enough to contain the rapidities of both observed gluons. 
The operators defined this way satisfy the standard com¬ 
mutation relations in the transverse momentum space: 

K(k),al j (p)] = (2tt fS ab 6^S^(k-p). (9) 

The integral over momenta k \, k 2 in eq. (|5]) contains a 
contribution from the region k\ = k 2 . If the wave func¬ 
tion |in) exhibits Bose enhancement, there is enhanced 
probability that the two gluons have equal momenta. 
This excess in the initial state will then translate into fi¬ 
nal state correlations. Note that this effect is suppressed 
by the squared number of colors 1 / N 2 , since Bose en¬ 
hancement is only operational for bosons with identical 
quantum numbers. 

We thus have to understand what is the nature of the 
projectile state | in), and in particular we need to calcu¬ 
late 

D(k 1 ,k 2 ) = {in\al\ki)al 3 (k 2 )a^(ki)a[(k 2 )\in). (10) 

Averaging over the projectile state in the standard 
CGC approach involves two elements. One needs to 
calculate the average over the soft degrees of freedom, 
as well as that over the valence color charge density. 
Conventionally this is done in the spirit of the Born- 
Oppenheimer approximation, namely first one averages 
over the soft gluon degrees of freedom at fixed valence 
color charge density p, and subsequently averages over 
the valence density distribution. 


The wave function of the soft fields for fixed valence 
color charge density for a dilute projectile is a simple 
coherent state 

\in) p = exp [<4 l (fc) + <(-&)]} 1°), ( n ) 

with the Weizsacker-Williams field b l a (k) = gpa^k)^. 

The averaging over the soft degrees of freedom leads 
to the well known expression for the observable in terms 
of the charge density: 

D{k u k 2 ) p = &i(* 1 )^(fc)&S(-fci)&£(-fc 2 ). (12) 

Since at fixed p, the soft gluon state is a coherent state, 
this expression does not seem to exhibit Bose enhance¬ 
ment. This is however misleading, since averaging over 
p is part of the quantum averaging over the initial state 
wave function | in). It is therefore instructive to reverse 
the conventional order of averaging, and average over the 
valence degrees of freedom first. The result of such a pro¬ 
cedure is a density matrix on the soft gluon Hilbert space. 
The subsequent averaging over this density matrix is a di¬ 
rect way to find out whether the projectile wave function 
exhibits Bose enhancement. 


C. The soft gluon density matrix 

The soft gluon density matrix of course depends on 
the weight for the valence color charge density. For il¬ 
lustrative purposes we choose the same Gaussian weight 
used in the glasma graph calculation, the McLerran- 
Venugopalan model |l5| , 

<•■■)„ =U j D[p] ... e~ fk ^) Pa(k)Pa( ~ k) , (13) 

where M is the normalization factor. 

Thus the density matrix of the soft gluons is given by 

p=M I D[p\e~ fk (14) 

x e */, |o)(0| e -i ^ 
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where we have defined <j> % a {k ) = a l a (k) + a£*(— k). The 
integral over p can be performed with the result 


p = e ~S k s ~t^ lkikj ^(- fe ) 


(15) 
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The interesting correlator is given by 
D(k 1 ,k 2 ) = tr\pal l (ki)al 3 (k 2 )a^(k 1 )a l b (k 2 )]- (16) 


It is a matter of straightforward algebra to show that 


tr[pal l (k)a J b (p)\ = (2tt ) 2 6 ab S { 2 ) (k-p) g 2 p 2 {p) ^ , 

tr[pa l a {k)a J b (p)\ = tr[pa}l (k)a \ 3 {p)\ 

= — (2n) 2 6 ab fi( 2 \k +p) g 2 p 2 {p) (17) 

P 4 

and then find 


D(k l ,k 2 ) = S 2 {N 2 - 1) 


2^r2 ^ 2 k l k l k 2 k 2 g 4 p 2 (kl)k- 2 (k 2 ) 


1 + 


1 


S(N 2 - 1) L 


U2 U2 U2 I„2 

rv-j^ '*'2 A/-j^ a^2 

S { 2 \k! - k 2 ) + 6 {2 \h + k 2 ) 


(18) 


In order to get eq. m, we have made substitutions of 
the type (2tt ) 2 S ( ' 2 \ki — k\) —> S, where S is the trans¬ 
verse area of the projectile. This regularization amounts 
to taking into account the discreteness of the transverse 
momentum spectrum of confined gluons. 


III. DISCUSSION 

The first term in eq. (1181) is the ” classical” term equal 
to the square of the number of particles. The second term 
is the typical Bose enhancement term, suppressed with 
respect to the first ” ‘classical” ’ term by the total number 
of degrees of freedom (color and area). The third term is 
specific to the density matrix at hand and, as explained 
in [7j, appears due to reality of the gluon field scattering 
amplitude. This establishes our point that the soft glue 
density matrix exhibits Bose enhancement, so that the 
likelihood of finding two gluons with the same transverse 
momentum is higher than average. Note that this effect 
is naturally subleading in N c as the enhancement is only 
effective if both gluons are in the same color state. 

As a typical Bose enhancement contribution, the sec¬ 
ond term in eq. (1181) is nonvanishing only when the mo¬ 
menta of the two gluons are equal. Note however that 
jfei and k 2 are not the momenta of observed gluons, but 
rather the momenta of gluons in the wave function of 


the incoming projectile. The two gluons then scatter 
on the target and acquire momenta p and q with the 
probability N(j> — k\)N(q — k 2 ), as indicated in eq. ©• 
Nevertheless it is clear that in favorable kinematics the 
initial state correlation of eq. m also appears as a cor¬ 
relation between the final state particles. Consider for 
example a situation where the incoming projectile wave 
function has an intrinsic saturation momentum Q s , and 
the momenta p and q are chosen to be of the same or¬ 
der as Q s i.e. p 2 ~ q 2 ~ Q 2 . In such a situation the 
production amplitude is dominated by the contribution 
from |fci|, \ k 2 \ ~ Q s , since the gluon density of the pro¬ 
jectile is dominated by those gluons close to saturation 
momentum. The delta function in the Bose enhancement 
term in eq. m is then smeared when convoluted with 
the scattering probability N(p — ki)N(q — k 2 ), but posi¬ 
tive angular correlations between the directions of p and 
q clearly survive. Thus the initial correlation is transmit¬ 
ted to the final state, provided fragmentation and final 
state effects in p-p and p -A collisions are small. On the 
other hand, when \p\, |qj Q s , the initial correlation is 
smeared out by the large momentum transfer from the 
target, and the correlation in the final state should dis¬ 
appear. These qualitative features are of course borne 
out by the numerical calculations of Cl¬ 
one interesting question naturally follows on from the 
above discussion. Fermions in the initial state wave func¬ 
tion surely experience Pauli blocking. One therefore may 
expect negative correlation between final state hadrons 
that originate from quarks or antiquarks in the initial 
state. Such correlation should exhibit anticollimation 
rather than collimation, and therefore a valley rather 
than a ridge at A 4> = 0. Whether such a valley extends 
to large relative rapidities between the observed particles 
is a question that should be explored. Quark-antiquark 
pairs are present in the hadronic wave function within 
the CGC approach at the next to leading order in a s via 
splitting of gluons. Since the gluonic wave function is 
boost invariant, the same is true for the quark and anti¬ 
quark distribution. However, the main question here is 
whether the fluctuations around some ’’mean field” are 
not too large to mask the correlations in rapidity event by 
event. Another way of saying it, is to recall that in our 
discussion of gluons only a single rapidity independent 
mode of the quantum gluon field was large in the CGC 
wave function. As a result any correlation extended over 
large intervals in rapidity. Whether a similar effect dom¬ 
inates the quark wave function has to be investigated. 
Work on these questions is ongoing [l6| . 

Perhaps an even more pressing question to understand 
is whether such valleys can be observed experimentally, 
given that the quark contribution is suppressed by a s 
relative to that of gluons. Here we see two possible av¬ 
enues. One point is that, as opposed to gluon contribu¬ 
tion to correlations, the quark contribution is not sym¬ 
metric under A <fi —> —A(f>. It thus can generate a non¬ 
vanishing V 3 coefficient within the CGC approach. Such 
mechanism will be quite different from the one explored 
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in 0 based on the idea of local anisotropy suggested in 
0. Another possibility is to trigger on final states which 
predominantly arise from quarks. For example it may be 
interesting to study correlations between two D-mesons 
(or B-mesons), since open charm (or beauty) should have 
a relatively larger component coming from fragmentation 
of quarks, rather than that of gluons 13- 
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APPENDIX 


diagrams 8a, 8 b and 8c were included in the numerical 
calculations of 13 for completeness, they are subleading 
at high momentum and do not contribute significantly 
to ridge correlations. We will therefore disregard these 
contributions in the following. 

Recall that we are calculating 


{out\al l {p)al J (k)a l a (p)a J b (k)\out) = J DpW[p\ 

. (19) 


We now transform the creation and annihilation op¬ 
erators by the sequence of unitary transformations 
e ~ l fk b a( k )[ a t The first transformation 

shifts the creation and annihilation operators by the 
Weizsacker-Williams field b l a {k), while the second one 
color rotates the gluon operators as well as the color 
charge density operators locally in the transverse space. 
In coordinate space representation, 

gtg*/* K( k )[ a l'( k )+ a a(- k )\ a 1 k h # fe )[ a “”( fe ) +a # -fc )],S' 
gi 

= S ab (x)a f b l (x) - - y)S ab (y)p b {y ). (20) 

We now use the fact that when acting on the state | in) p , 
the charge density operator and the gluon creation (and 
annihilation) operators are identical in the sense that 
\d l a^{x) — p a {x)\ | in) p = 0. We can thus write in mo¬ 
mentum space 


gtgi/fc Ki k )[ af a( k )+ a a(- k )] a ^ i ^p'j e - i Ik b 'a( k )[ a i'( k )+ a a(- k )]cj 


a b {ki)\in) p . (21) 


/ S ab (p-ki) 

i 

1 

i_ 

) k\ 

p 


Here we explain in more detail the meaning of fig. 
[T| The diagrams are depicted before averaging over the 
projectile wave function. It is however straightforward 
to identify the correspondence between the diagrams on 
fig. IH and the glasma graphs of [3- The two possi¬ 
ble nontrivial Wick contractions between a’s and a'’s 
in the diagram of Type A generate the graphs la and 
7b of [3- The contractions (a t (fci)a(fc 3 ))(a' l '(fc 2 )a(fc 4 )) 
in Type B and C generate the glasma graphs (not 
explicitly drawn in [13() which have the target and 
projectile charge densities interchanged. They can 
be interpreted as Type A but with gluons originating 
from the target rather than from the projectile. The 
contraction (a^(k±)a(k 4 ))(a^(k 2 )a(k 3 )) in Type B and 
(aHki )aUk 2 )){a(k 3 )a(k 4 .)) in Type C diagrams generate 
diagram 8c and 8 b of 1131 ] respectively, whose contribution 
is proportional to S^ 2 \p q= q). As discussed in 0 , these 
diagrams are of the type that lead to HBT correlations if 
the translational invariance approximation is relaxed. Fi¬ 
nally the contraction (a* (ki)a' (k 2 ))(a(k 3 )a(k4,)) in Type 
B and (a* (fei )a(k/i))(a ' fejaffca)) in Type C give the dia¬ 
gram 8 a of jl.3j . The contribution of this diagram, as dis¬ 
cussed in [13[ , is suppressed at high momenta. Although 


Introducing the notation 


Sab(P’ k ) 


Sab{p - k ) 



( 22 ) 


we can now write the two gluon inclusive cross section 
in terms of the expectation value in the initial projectile 
state as 

[ (S^ (p, h)SC 2 (q, k 2 )S2 3 (-p, -k 3 )S* (- q , -fc 4 )> 

x (HOoi 1 ( fc i) a a 2 2 (k 2 )a x a 3 3 (k 3 )a l a \(k 4 )\in). (23) 

The next step is to average over the target fields. Within 
the glasma graph model one assumes that the averages 
of the S-matrices are translationally invariant, pairwise 
factorize, and the color structure of the averages is the 
same as in the dilute limit: 


<■ Sab{p)S cd (k )) = r be f cde N(p)5(p + k ). (24) 

In this approximation the averaging over the target fields 
can be performed. This averaging leads to graphs of three 
distinct topologies depicted in fig. |I] The blobs denote 
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the appropriate vertices, as in eq. ([22]). For example, the two distinct blobs in the Type A diagram are S 

and S ij - 
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